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This paper presents the modeling and simulation for a prototype of the solar-powered Stirling engine 
working at the low temperature range. A mathematical model for the thermal analysis of the solar- 
powered low temperature Stirling engine with heat transfer is developed using Matlab program. The 
model takes into consideration the effect of the absorber temperature on the thermal analysis like as radi¬ 
ation and convection heat transfer between the absorber and the working fluid as well as radiation and 
convection heat transfer between the lower temperature plate and the working fluid. Hence, the present 
analysis provides a theoretical guidance for designing and operating of the solar-powered low tempera¬ 
ture Stirling engine system, as well as estimating output power from the solar Stirling engine using 
Matlab program. This study attempts to demonstrate the potential of the low temperature Stirling engine 
as an option for the prime movers for Photovoltaic tracking systems. The heat source temperature is 
40-60 °C as the temperature available from the sun directly. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Solar energy is the energy of sunlight collected and used to pro¬ 
vide electricity, heating, cooling and other purposes. This clean en¬ 
ergy reduces the consumption of fossil fuels and also reduces the 
environmental pollutants [1 . Solar energy is one of the more 
attractive renewable energy sources that can be used as an input 
energy source for heat engines. In fact, any heat energy source 
can be used with the Stirling engine [2 . 

The direct conversion of solar power into mechanical power re¬ 
duces both the cost and complexity of the prime mover. Stirling 
engines are able to use solar energy that is a cheap source of en¬ 
ergy. The increasing interest in Stirling engines is largely due to 
the fact the engine is more environmentally friendly than the 
widely used internal combustion engine [3 . Stirling engines are 
thermodynamic devices working theoretically on the Stirling cycle 
and using a compressible fluid, such as air, hydrogen, helium, 
nitrogen or even vapors, are used as working fluids. The Stirling en¬ 
gine offers possibility for having high efficiency engine with less 
exhaust emissions in comparison with the internal combustion 
engine. The earlier Stirling engines were huge and inefficient. In 
recent years, a number of new Stirling engine models have been 
developed to improve the deficiencies [2 . 
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Tlili et al. [4] presented the engine works with a closed cycle 
and used several gases as working fluid. Glushenkov et al. [5] eval¬ 
uated thermodynamic analysis of an unconventional heat engine. 
The engine studied has a number of advantages compared to state 
of the art Stirling engines. Cheng et al. [6] estimated theoretical 
analysis of the effects of the geometrical parameters on the shaft 
work of the Stirling engines. Campos et al. [7] showed a mathemat¬ 
ical model to simulate the operation of Stirling engines. A thermo¬ 
dynamic optimization of the cycle for maximum thermal efficiency 
was conducted. Iskander Tlili et al. [8] estimated maximum power 
and efficiency at the maximum power point of an end reversible 
Stirling heat engine with finite heat capacitance rate of external 
fluids in the heat source/sink reservoirs with regenerative losses 
are treated. Cheng et al. [9] presented development of a dynamic 
model for the beta-type Stirling engines with rhombic-drive mech¬ 
anism. Cheng et al. [10] presented numerical predictions of rela¬ 
tionship between operating speed and shaft power output of 
Stirling engines. 

Aksoy and Cinar [11] investigated the effect of the grooved 
displacer cylinder on engine performance from kinematic and ther¬ 
modynamic points of view. Chen et al. [12] studied experimentally 
a helium charge y-type twin power piston Stirling engine to under¬ 
stand the effects of several regenerator parameters on the overall 
performance of the engine. Formosa and Despesse [13] presented 
an analytical model for design and optimization of FPSE. To deal 
with the discrepancy between the high theoretical efficiency of 
Stirling and a constructed prototypes the major losses are be taken 
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Nomenclature 




SL 

standard longitude, min 

Fa } - Pl 

the view factor coefficient between point al & 

LL 

local longitude, min 


Plate 1 

L 

local latitude, min 

Fp 2 a 2 

the view factor coefficient between Plate 2 & 

6 

the declination angle, ° 


point a2 

h 

the hour angle, ° 

Apt 

the Area of the Plate 1, m 2 

0) 

the Zenith angle, ° 


the Area of the Plate 2, m 2 

0 

the solar incidence angle, ° 

G 

the Stefan-Boltzmann constant (5.669 x 

p 

the surface tilt angle from the horizontal, ° 


10 8 W/m 2 K 4 ), W/m 2 K 4 

Ct 

the global irradiance, W/m 2 

£ Pi 

the Plate 1 emissivity 

^B t 

total beam solar radiation on a tilted surface, 

£ Pi 

the Plate 2 emissivity 


W/m 2 

h c 

the connective heat transfer coefficient, W/m 2 I< 

^Dt 

total diffuse solar radiation on a tilted surface, 

K 

the thermal Conductivity of Fluid (air), W/m K 


W/m 2 

L 

the length of fluid layer = 1.5 cm = 0.015 m, m 

GR t 

total ground-reflected solar radiation on a tilted 

N„ 

the Nusselt Number 


surface, W/m 2 

G r 

the Grasshof number 

C B 

beam radiation on a horizontal surface, W/m 2 

Pr 

the Prandtl number 

Cb„ 

beam radiation in the direction of the rays, W/ 

a, n 

the constants based on fluid circulation 


2 

m 

P T 

the coefficient of thermal expansion, 1 /K 

Rb 

the beam radiation tilt factor 

P 

the dynamic viscosity of Fluid (air), kg/m s 

Rd 

the diffuse radiation tilt factor 

P 

the density of Fluid (air), kg/m 3 

Rr 

the ground-reflected radiation tilt factor 

g 

the gravity of acceleration = 9.8 m/s 2 , m/s 2 

p 

the ground reflectance, % 

(T,-T 2 ) 

the temperature difference, K 

Tair 

the temperature of the container air of Stirling 

C P 

the specific heat, kj/kg K 


engine, I< 

r 0 Xi 

outer radius & inner radius of Plate 1, m 

T Pl Pi 

the temperature of upper and lower layer of 

L 

length of air fluid inside heating flow, m 


Stirling engine, I< 

Di, D 2 

diameter of Plate 1 & Plate 2, m 

Ta, 

the temperature of the surrounding air in sun 

K ■ 

MX air 

the conductivity of air, W/m K 


area above Plate 1, I< 

Pair 

the viscosity of air, kg/m I< 

Ta 2 

the temperature of the surrounding air below 

P 

the power output of the Stirling engine, watts 


Plate 2, K 

Pm 

the mean average gas pressure in the container 

d 

u Pi Pi air 

the thickness of each layer, m 


of Stirling engine, (Pa) or (MPa, if volume is in 

Pp\ ,p 2 ,air 

the density of each layer, kg/m 3 


cm 3 ) 

Cp p u p 2 ,air 

the heat capacity of each layer, J/kg K 

V 

the swept volume of the expansion space in the 

^Pi Pi 

the amount of solar radiation falling on the sur- 


container of Stirling engine, (m 3 ) or (cm 3 , if 


face of each layer, W/m 2 


pressure is in MPa) 

c t 

the total solar radiation on a surface, W/m 2 

F 

the engine cycle frequency of Stirling engine 

a Pl 

the absorptance of Plate (1) = 0.95 (Plate Black 


container, Hz 


Paint) 

Bn 

the Beale number depend on air heater temper- 

a P2 

the absorptance of Plate (2) = 0.25 (Plate Silver 


ature of Stirling engine 


Paint) 

Pair 

the density of air (=1.127 kg/m 3 ), kg/m 3 

h C ^-p^h C Arp 2 -a 2 

the convection heat transfer coefficients out- 

R 

the constant of proportionality, is known as the 


side Stirling, W/m 2 I< 


universal gas constant (=287 J/kg I<), J/kg I< 

^c,lr Pl -air? hc,air-p 2 

the convection heat transfer coefficients inside 

Tair 

the temperature of the container air of Stirling 


Stirling, W/m 2 K 


engine, I< 

^Pi Pi 

the view factor coefficient between Plate 1 & 




Plate 2 




into account. Ahmadi et al. [14 presented a thermodynamic model 
to determine output power and thermal efficiency of the solar-dish 
Stirling heat engine. A final optimal solution was selected from 
available solutions of the Pareto frontier using three decision mak¬ 
ing methods including LINMAP, TOPSIS and fuzzy methods. Then, 
Ahmadi et al. [15] presented finite-time thermodynamics to deter¬ 
mine output power, thermal efficiency and rate of entropy gener¬ 
ation of the Solar Dish-Stirling engine and factors such as dish 
collector performance, conductive and radiative heat transfer 
mechanism at heat source and heat sink of the engine are involved 
in the analysis. 

Sripakagorn et al. [16] developed a prototype Stirling engine 
working at the moderate temperature range. Abdullah et al. [17] 
presented the design considerations to be taken in designing of a 
low temperature differential double-acting Stirling engine for solar 
application. Li et al. [18] developed a mathematical model for the 


overall thermal efficiency of the solar-powered high temperature 
differential dish-Stirling engine with finite-rate heat transfer, 
regenerative heat losses, conductive thermal bridging losses and fi¬ 
nite regeneration processes time. Kongtragool et al. [19] provide a 
theoretical investigation on the thermodynamic analysis of a Stir¬ 
ling engine. Kongtragool et al. [2] provided a literature review on 
solar-powered Stirling engines and low temperature differential 
Stirling engines technology. 


2. Design of solar-powered low temperature Stirling engine 

In this paper, a low temperature Stirling engine has been devel¬ 
oped. This low temperature Solar-powered Stirling motor design 
was conceived for the transduction of solar energy at a necessary 
low temperature difference of 10:20 °C as shown in Fig. 1. During 
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Fig. 1. Stirling engine 3-D design. 



Fig. 2a. Stirling engine (Heating Cylinder) design. 



Fig. 2b. Stirling engine (Heating Cylinder components) design. 1. Elbow Joint, 2. 
Power Cylinder, 3. Top Plate, 4. Heating Cylinder, 5. Lower Plate. 


/V 

2.5 cm 


5.5 cm 


v v 


T 

4 cm 

.1 


1.5 cm 


h-H 

11 cm 

Fig. 2c. Stirling engine (Heating Cylinder dimensions) design. 



the construction of the low temperature solar-powered Stirling en¬ 
gine, an emphasis was put on temperature differences between 
two plates upper and lower, which is the important condition for 
operating the engine. 

This configuration uses separated cylinders for the displacer 
and the power pistons, with the power cylinder connected to the 
displacer cylinder. The displacer moves working fluid between 
the hot space and the cold space of the displacer cylinder (Heating 
Cylinder). In this configuration, the power piston both compresses 
and expands the working fluid [2,20]. 

A solar-powered Stirling engine has a single power piston (26) 
as shown in Figs. 3a-3c arranged within the power cylinder (2) 
in Figs. 2a-2c on the same shaft as a displacer piston (27) in 
Figs. 3a-3c. The displacer piston is a loose fit and does not extract 
any power from the expanding gas but only serves to shuttle the 
working gas between the hot and cold heat exchangers. When 
the working gas is pushed to the hot end of the cylinder it expands 
and pushes the power piston. When it is pushed to the cold end of 
the cylinder it contracts and the momentum of the machine, usu¬ 
ally enhanced by a flywheel, pushes the power piston the other 
way to compress the gas. Unlike the alpha type, the beta type 
avoids the technical problems of hot moving seals. 

In this way, the movement of the working piston, discharger 
and crank mechanism can be visualized. The sealing of the heat 
transduction space is ensured by two O-rings, which simulta¬ 
neously serve to secure the upper and lower plates with the 


casing wall. Assembly is achieved by a re-sealable Snap-on 
connection. 

To minimize friction, the crankshaft (13) as shown in Figs. 3a- 
3c and both connection rods (3, 4) as shown in Figs. 4a and 4b 
are mounted in precision miniature ball bearings (23). Since the 
working piston runs in a dry cylinder, the low temperature solar- 
powered Stirling engine is extremely maintenance resistant. 

The upper and lower Plates (3 and 5) were made from 0.3 cm 
thick aluminum plate. The plates were then machined on a Bridge¬ 
port mill, using a rotary table. To maintain concentricity, the upper 
and lower plates were machined simultaneously as shown in 
Figs. 2a-2c. 

The displacer (27) is a special-purpose piston, used in this type 
of solar Stirling engine, to move the working gas back and forth be¬ 
tween the hot and cold heat plates. Depending on the type of en¬ 
gine design, the displacer is sealed to the cylinder; it may be a 
loose fit within the cylinder, allowing the working gas to pass 
around it as it moves to occupy the part of the cylinder inside. 
The displacer (27) itself was made from 0.6 cm thick high density 
foam board, much like the blue foam used for insulation in the con¬ 
struction of houses. The foam disk was easily formed into place 
with the help of a belt sander. Proper thickness of 0.8 cm was ob¬ 
tained by a band saw, and cleaned off by the belt sander as shown 
in Figs. 3a-3c. 

The displacer rod (21) was supposed to be made out of steel, but 
due to limited resources the only 0.127 cm rod that was available 
was brass-welding rod as shown in Figs. 3a-3c. 
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Fig. 3a. Stirling engine (rotating part) design. 



Fig. 3b. Stirling engine (rotating part components) design. 1. Flywheel, 2. Crank¬ 
shaft, 3. Shim, 4. Ball Bearing, 5. Bearing Bush, 6. Ball Bearing, 7. Shim, 8. Crank Disc 
1, 9. Socket, 10. Shim, 11. Crank Pin, 12. Shim, 13. Crank Disc, 14. Shim, 15. Ball 
Bearing, 16. Shim, 17. Head Screw, 18. Connecting Rod, Short, 19. Adjusting Collar, 
20. Connection Cylinder, 21. Displacer Rod, 22. Displacer Bushing, 23. Ball Bearing, 
24. Connecting Rod, Long, 25. Cylinder Pin, 26. Power Piston, 27. Displacer. 


The bearing plate (5) was cut out of the same piece of aluminum 
that the top and bottom plates came from. The bearing housing (4, 
6) and collar (19) were made from aluminum using a lathe. Other 
parts such as the crankpin (11), flywheel hub (1), crankshaft (8, 
13), and displacer crankpin (15), were also turned on a lathe as 
shown in Figs. 3a-3c. 

One of the most difficult pieces to machine was the Displacer 
Chamber Ring. It was a rough-cut Plexiglas ring with an inside 
diameter of 10.5 cm. The height of the ring had to be between 
1.4 cm and 1.5 cm. An improper height would offset the volume in¬ 
side the chamber. To aid in the making this piece, a template was 
made out of plywood. Using a Rotozip tool with a circular attach¬ 
ment, a path 0.5 cm was cut. The ring was then placed into the 



k- 

9 cm 


Fig. 3c. Stirling engine (rotating part dimensions) design. 



Fig. 4a. Stirling engine (connection part) design. 



^II 

a j 


Fig. 4b. Stirling engine (connection part components) design. 1. Connection Ball 
Bearing and Rings, 2. Head Screw, 3. Connecting Rod, Short, 4. Connecting Rod, Long. 


groove and sanded smooth on a belt sander. As the plastic was 
sanded, the excess material seemed to melt away. The ring had 
to be perfectly flat on the top and bottom faces, to keep air from 
leaking around the ring when it was sandwiched between the 
upper and lower plates as shown in Figs. 2a-2c. 

The flywheel (1) was made by adhering three compact disks to¬ 
gether with epoxy. The rods connecting the piston (26) to the right 
bearing assembly (23), and the displacer rod (21) to the left bearing 
assembly. This material is rigid enough that it will not deform un¬ 
der normal usage, but easy enough to shape into the proper form 
using simple tools as shown in Figs. 3a-3c and 4a-b. 
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Fig. 5. Beam radiations on horizontal and tilted surfaces. 


Table 1 

Coefficients for 1990 Perez diffuse radiation model [23 . 


£ 

/n 

f\2 

/l3 

h\ 

f*22 

/23 

1-1.1065 

-0.008 

0.588 

-0.062 

-0.060 

0.072 

-0.022 

1.065-1.230 

0.130 

0.683 

-0.151 

-0.019 

0.066 

-0.029 

1.230-1.500 

0.330 

0.487 

-0.221 

0.055 

-0.064 

-0.026 

1.500-1.950 

0.568 

0.187 

-0.295 

0.109 

-0.152 

0.014 

1.950-2.800 

0.873 

-0.392 

-0.362 

0.226 

-0.462 

0.001 

2.800-4.500 

1.132 

-1.237 

-0.412 

0.288 

-0.823 

0.056 

4.500-6.200 

1.060 

-1.600 

-0.359 

0.264 

-1.127 

0.131 

>6.200 

0.678 

-0.327 

-0.250 

0.156 

-1.377 

0.251 





Ta* 

Fig. 6. Simplified system diagram with Stirling engine heating container and 
temperature in sides. 


Table 2 

Estimation of a and n constants. 


Surface 

(Gr-Pr) 

a 

n 

Vertical plates/cylinders 

10 4 to 10 9 

0.59 

0.25 


10 9 to 10 12 

0.13 

0.33 

Horizontal pipes 

10 3 to 10 9 

0.53 

0.25 

Horizontal plates 

10 5 to 2 x 10 7 

0.54 

0.25 

Heated face up or cooled face down 

2 x 10 7 to 3 x 10 10 

0.14 

0.33 

Vertical plates heated face up or cooled 

3 x 10 5 to 3 x 10 10 

0.27 

0.25 

face down 






Heftier temperature (K) 


Fig. 7. Beale number as a function of source temperature [26 . 


algebraic sum of three solar components; first is the direct (beam) 
solar radiation, the second is the diffuse solar radiation, and the 
third is the ground reflected solar radiation [21 . A flat surface ab¬ 
sorbs beam (G Bt ) diffuse (G Dt ) and ground-reflected (G Rt ) solar radi¬ 
ation; that is [22 , 

Q = f^Bt + Got + Gr,. ( 1) 


There are two positions in move Stirling engine: 


3.1. Beam radiation 


(1) When the displacer is near the bottom of the Heating Cylin¬ 
der; inside the engine most of the gas has been heated by the 
heat source (sun energy) and it expands. This causes the 
pressure to increase which forces the piston down. 

(2) When the displacer is near the top of the Heating Cylinder; 
most of the gas in the engine has now cooled and contracts 
causing the pressure to decrease, which in turn allows the 
piston to move up and compress the gas. 

3. Solar radiation modeling on tilted surfaces calculations 

In this section a solar simulation has been made to estimate 
average global radiation using hours of bright sunshine and calcu¬ 
lated solar radiation data available for solar powered Stirling en¬ 
gine in Egypt. In this part, an empirical model is proposed for 
estimating daily global solar radiation on a horizontal and on a 
tilted, south-facing surface by the day of the year. Also, presents 
the implementation of a generalized solar radiation simulation 
model using MATLAB/GUI interface. Global solar radiation is the 


The beam component is calculated using the direct normal irra- 
diance (G Bn ) from the simulation model. 

As shown in Fig. 5, the beam radiation on a tilted surface is 

Gb, = C B „ cos(0) (2) 

and on a horizontal surface, 

G b = G Bn cos(<D) (3) 

It follows that 


_ c Bt _ cos(e) 

G B ~ cos(O) 

The beam radiation component for any surface is 


G Bt 


G r R 


B^B 


( 4 ) 

( 5 ) 


In Eq. (4), the zenith angle (O and the incident angle (0) can be 
Simplified. Therefore, Eq. (4) becomes 

_ cos(0) _ sin(L - p) sin(S) + cos(L - p) cos(S) cos(h) 

B ~ cos(O) - sin(L) sin(S) + cos(L) cos(5) cos(h) 
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Solar Radiation at Stirling Engine 



Time (hrs) 

Solar Radiation at Stirling Engine 


Fig. 8. Practical solar radiation data in Stirling engine heating container during day hours on 2 March 2013 in Giza. 


Table 3 

Data needed for solar radiation calculations in Giza city for Matlab simulation model 

[27]. 


Location 

Egypt - Giza National Research Centre 

The standard meridian, Egypt 

31°17' East 

The local longitude, Giza 

31°l(y East 

The local latitude, Giza 

TOW North 

Date 

Changes over the year 

Time 

Changes over the day 

The ground reflectance 

0.27 

Tilt angle, p 

30 


The ratio of beam radiation on the surface to that on a horizon¬ 
tal surface at any time is given by 


\](\ - cos 2 (8) sin 2 (h)) 
sin(L) sin(S) + cos(L) cos(S) cos(h) 



3.2. Diffuse radiation 


The diffuse radiation on a tilted surface is: 


Go t — Gd^d 
Gp t = (0.11xG Bn )i?D 


( 8 ) 

(9) 


Three reflected radiation models were implemented previously: 
the first generation model of Liu and Jordan, the second generation 
irradiance model of Klucher, and the third generation 1990 model 
of Perez et al. (see Fable 1 [23 ). The three models differ in the com¬ 
putation of R d This Matlab software using first model (Isotropic 
Model). 

In this model, the calculations of R D is simple and given by 



1 + cos(p) 
2 



The diffuse radiation component of the Liu and Jordan [11,20] 
isotropic model G Dt is 




|U+_cosym 



3.3. Ground reflected radiation 

The ground reflected radiation on a tilted surface is: 

G Rt = (G b + G d )R r (12) 

The solar radiation simulation model and results will be pre¬ 
sented in Section 6. 


4. Thermal analysis modeling for solar-powered Stirling engine 


The Stirling hot air engine is a simple type of engine that uses a 
compressible fluid as the working fluid. Because the working fluid 
is in a closed system, there are no problems with contamination 
and working fluid costs. 

This part of paper presents the modeling and simulation for a 
prototype of the solar-powered Stirling engine working at the 
low temperature range. A mathematical model for the thermal 
analysis of the solar-powered low temperature Stirling engine with 
heat transfer is developed using Matlab program. The model takes 
into consideration the effect of the absorber temperature on the 
thermal analysis like as radiation and convection heat transfer be¬ 
tween the absorber and the working fluid as well as radiation and 
convection heat transfer between the lower temperature plate and 
the working fluid as shown in Fig. 6. Hence, the present analysis 
provides a theoretical guidance for designing and operating of 
the solar-powered low temperature Stirling engine system, as well 
as estimating output power from the solar Stirling engine using 
Matlab Program as shown in Fig. 19. This program attempts to 
demonstrate the potential of the low temperature Stirling engine 
as an option for the prime movers for Photovoltaic tracking sys¬ 
tems. The heat source temperature is 40-60 °C as the temperature 
available from the sun directly. The details of these results are pre¬ 
sented in next section. 

Eqs. (13)-(15) are a general relationship for energy balance 
equations for three points at plate 1 and plate 2 and at the con¬ 
tainer midpoint [24 . 


a Pl G t + h ca 1 - Pl (T fll — T Pl ) + {7£ Pl i 7 a 1 _ Pl (T fli - T p2 ) 

= hc,lrpi-air(Tp^ ~ T a ir) + ^Pi^Pi p 2 ~ ^p 2 ) 
h C ,air-p 2 (Tair ~ Fp 2 ) + <7£ Pl F Pl _ P2 (A P| / Ap 2 ) (T pi — T pi ) 

= hc,lrp 2 -a 2 (T P2 — Ta 2 ) + (J£p 2 F p2 _ a2 (T p 4 — T a 4 ) 


h 


c,p^-air\ M pi 


(Tp T air ) — h cair p2 (T air T p ) 


(13) 

(14) 

(15) 
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The Stefan-Boltzmann constant o = 5.669 x 10~ 8 W/m 2 K 4 . The 
specific heat of plate 1 and plate 2 Cp p i, Cp p2 such as for aluminum 
frame is 897 J/kg K. The specific heat Cp air such as for air inside cyl¬ 
inder is 1005J/kgK. The solar absorptance coefficients for surface 
materials plates a Pl , a P2 are 0.95 and 0.25. The density for alumi¬ 
num p pl and p p2 is 2712 kg/m 3 and for air p air are 1.205 kg/m 3 . 
The thickness of each layer, d Pl P2 air are 0.001, 0.001, and 0.015 m. 

Convective heat transfer is one of the major modes of heat 
transfer and convection is also a major mode of mass transfer in 
fluids. The general equations for thermal coefficients as connective 
heat transfer, Nusselt number, Grasshof number, Prandtl number 
coefficients are shown in Eqs. (16)-(19) [24]. 


h c 

N h 

Gr 

Pr 


= a(G r xP r ) n 
gPp 2 L 3 (T i 

p 2 


Cp/I 

K 



( 16 ) 

(17) 

(18) 
(19) 


The general equation for Nusselt number equation at natural 
convection is shown in Eq (17). Convection heat transfer values 
are very specific to the geometry of the surface and the heat 


a sr 



SoCar (Radiation (Program 



i 


Fayoum 'University 

Measuring Pfa.ce 


Standard longitude 

31 Degree ~~\J I Minute q Second 

Locaf longitude 

31 'Degree 10 Minute q Second 

Locaf Latitude 

30 Degree Q Minute q I Second 


r Measuring Date 

Day MontR Year 

2 3 2013 


Sfationaf Research Centre 

r~ Figure addition - 


Totalsofar radiation u, it H Day time 


i"Required Input- 


r Measuring Pface Time - 

Locaf standard time 

12 Hour ~o I Minute 0 Second 


(ground Type 

Cement ▼ 

THe Pface Direction 
from greniucR 

‘East ▼ 

TRe Lifted Angle 

30 



Output Data of So far (Radiation at standard Time- 


Output 


Lota f So far Radiation 
Seam solar radiation 


1012.52 <w/ / m2 

850.213 W7m2 

Diffuse solar radiation 140.115 tV/m2 

ground-reflectedsolar radiation 22 1899 W/m2 


Programmed Sy E^rg jlfimed Za£ana Hafez 


Output Data ofTotafSofar Radiation per day- 


Lotaf Sofar Radiation per day 

5.68633 XyWt/m2 per day 

Sunrise rime 6 Hour 18 Minute 27 Second 

Sunset Time 17 Hour 41 Minute 33 Second 


Contact Ids 


(a) 





Time (hrs) 




■■^•Total Solar Radiation 


•Beam Solar Radiation 


►Diffuse Solar Radiation 


Ground Reflected Solar Radiation 



Fig. 9. Estimation of daily solar radiation on an inclined plane (south facing - Tilted 30°) during day hours in Giza on 2 March 2013. (a) Matlab GUI/Interface simulation 
model (b) solar radiation data during day hours. 
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transfer conditions as shown in Table 2, we can estimate a, n con¬ 
stants values from Table 2 [24]: 

As shown in Fig. 2c for Stirling engine container plate 1, the out¬ 
er radius of plate 1 is 5.5 cm, and the inner radius hole of plate 1 is 
1 cm. The length (L) between plate 1 and plate 2 is 1.5 cm. The 
thickness of each layer, d Pl P2 air are 0.001, 0.001, and 0.015 m. 

The specific heat Cp air such as for air inside cylinder is 1.005 kj/ 
kg K. The density for air p a/r is 1.205 kg/m 3 . The solar absorptance 
coefficients for surface materials plates a Pl ,a P2 are 0.95 and 0.25. 
The thermal Conductivity of air K a ir is 0.0257 W/m K. The 


coefficient of thermal expansion of air p air is 3.43 1 /K. The dynamic 
viscosity of air /i air is 1.983 * 10~ 5 g/ms. The gravity of acceleration 
g is 9.8 m/s 2 . 

Eq. (16) is a general relationship for connective heat transfer 
coefficient. As shown in Eq. (20), it can be reduced to identified 
forms for four layers of Stirling engine container 

h ca , =0.895055 (T fl , -T Pl )°' 33 ,fi CJ , 2 -a 2 =0.903129(T P2 -T„ 2 f 33 , 

hcrt-air = 28.0346(T Pl -T alr ) n2 \h calr „ 2 = 28.0346 (T air - T Pl ) 025 

( 20 ) 
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Fig. 10. Estimation of daily solar radiation on a horizontal plane during year months at Giza, (a) Matlab GUI/Interface simulation model, (b) Solar radiation data during year 
months. 
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Coefficient s is the layer emissivity, and F is the view factor be¬ 
tween two surfaces. The emissivity of the upper aluminum plate is 
1. The emissivity of the lower plate for aluminum is 0.77. The solar 
absorptance coefficients for surface materials plates a Pl ,a P2 are 
0.95 and 0.25 [25]. 


X-Y 


P\ Pi 


0.894,F fllPl 


1 + cos(fi) _ 1 c _1 +cosQ?)_ 1 

— 1 5 r p 2 -fl 2 _ n — 1 


( 21 ) 


The view factor from the inner side for aluminum upper plate to 
the inner side for lower plate F Pl _ P2 is 0.894, and 1 from the outer 
surrounding air to the outer side for upper aluminum plate F ai _ Pl . 
The view factor from the outer side for aluminum lower plate to 
the outer surrounding air F Pl tt2 is 1 as shown in Eq. (21). 

As shown in Eqs. (13)-(15), the energy balance equations at 
plate 1 and plate 2 and at the container midpoint in the condition 
of steady state can be reduced too much simpler forms as Eqs. 
(22)-(24) as follow: 
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Fig. 11. Estimation of daily solar radiation on an inclined plane (South facing - Tilted 30°) during year months at Giza, (a) Matlab GUI/Interface simulation model, (b) Solar 
Radiation Data during year months. 
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0.95G t + 0.895055(T fll - T Pl )’ 33 + 5.669 x lO" 8 )!^ -T 4 f2 ) 

= 28.0346 (T Pl - T air ) h25 + 5.068 x lO' 8 )!^ - T 4 2 ) (22) 

28.0346(T afr - T P2 ) 125 + 4.9 x lO' 8 )^ - T 4 pj ) 

= 0.903129(T P2 - T„ 2 ) 133 + 4.36 x 10 _s (T£ 2 - T^) (23) 

Tair = Tpi 2 Tp2 (24) 

5. Sirling engine power output modeling using Beale formula 

In mechanical engineering, the Beale number is a parameter 
that characterizes the performance of Stirling engines. It is often 
used to estimate the power output of a Stirling engine design. 
For engines operating with temperature differential, typical values 


for the Beale number range from (0.11) to (0.15); where a larger 
number indicates higher performance. To estimate the power out¬ 
put of an engine, nominal values are estimated for the Beale num¬ 
ber (as shown in Fig. 7), pressure, swept volume and frequency, 
then the power is calculated as the product of these parameters, 
as follows: 

P = B n P m VF (25) 

Pm v = mRT air (26) 

P m = p air RT air = 323.49 T air (27) 

As shown in Fig. 2c for Stirling engine container plate 1, the out¬ 
er radius (r) of plate 1 is 5.5 cm, and the length (L) between plate 1 
and plate 2 is 1.5 cm. 

V = nr 2 L = 142.6 * lCT 6 m 3 (28) 

P = B„P m VF = 0.046129674 B„T air F (29) 


Stirling Engine Temperature Outputs 



Fig. 12. Estimation of Stirling engine temperature for (Plate 1, Plate 2, and Air in Stirling engine container) during day hours on 2 March 2013 at Giza. 


Stirling Engine Temperature Outputs 



Fig. 13. Stirling engine temperature for (Plate 1 and Plate 2) in Stirling engine container during day hours on 2 March 2013 at Giza. 
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6. Results 

Results from this study indicate that the engine performance in¬ 
creases with increasing solar intensity. The engine torque, shaft 
power, thermal efficiency, speed, and heating source temperature 
also increase with increasing solar intensity. In fact, it can be said 
that the maximum engine torque, shaft power, and thermal effi¬ 
ciency increase with increasing heating source temperature. The 
Beale number of this engine increases with decreasing tempera¬ 
ture ratio or with increasing heating source temperature. 

At the maximum solar intensity of 900 W/m 2 , or heat input of 
8.55 J/s, the engine produced a maximum shaft power of 
0.5171 W at 126 rpm, with a heating source temperature of 
64.2 °C and air temperature for heating container of 51.52 °C. 

Practical results and Simulation Matlab GUI/Interface model 
shows 

(1) Solar radiation data: solar radiation MATLAB /GUI interfac¬ 
ing Results. 


(2) Thermal analysis in low temperature solar-powered Stirling 
engine for temperature of two plates of container and air in 
container: 

i. Solar powered Stirling engine temperature results. 

ii. Matlab simulation program for solar powered Stirling 
engine temperature results. 

(3) Power outputs data for Stirling engine: 

i. Solar powered Stirling engine power estimation results. 

ii. Matlab simulation program for solar powered Stirling 
engine power estimation results. 


6.2. Solar radiation results 

Fig. 8 shows the amount of daily solar insolation during the 
daytime. It is shown in Fig. 8 that the maximum solar radiation ap¬ 
plied in Stirling engine is 900 W/m 2 at 12:00 PM. The results of dai¬ 
ly solar radiation during daytime are showed. 


Stirling Engine Temperature Outputs 



Time (hrs) 

Output Temperature in Plate 2 — Input Temperature in Plate 1 


Fig. 14. Input and output temperature for (Plate 1 and Plate 2) in Stirling engine container during day hours on 2 March 2013 at Giza. 


Stirling Engine Temperature Outputs 



Fig. 15. Plate 1 and input temperature for Stirling engine container during day hours on 2 March 2013 at Giza. 
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Input data for place site (Giza) as latitude and longitude are 
showed in Fable 3 [27] for Matlab simulation model. Fig. 9 shows 
Matlab/GUI estimation of Daily Solar Radiation on an inclined solar 
Stirling engine (South facing - Tilted 30°) during day hours on 2 
March 2013 at Giza, the results in Fig. 10 shows estimation of daily 
solar radiation on a horizontal plane during year months at Giza, 
and the results in Fig. 11 shows estimation of daily solar radiation 
on an inclined plane (South facing - Tilted 30°) during year months 
at Giza. 

6.2. Temperature results for Stirling engine 

The temperature measurement experiments took place on 2 
March 2013 from 07:00 AM to 5:00 PM by LM35 Sensor and data 


acquisition (NI USB 6008). During these experiments the weather 
conditions were very good and were almost no clouds in the sky 
else from 4:00 PM to 5:00 PM. The experimental results of Stirling 
Engine Temperature for (Plate 1, Plate 2, and Air in Stirling engine 
container) are shown in Fig. 12. Fig. 12 shows the temperature 
variations of each layer (upper and lower plates) during the day¬ 
time. The maximum temperature of air inside Stirling engine was 
51.52 °C at 1:00 PM; the results in Fig. 12 clearly show the effects 
of the upper and lower plates temperature for air temperature in¬ 
side heating container. 

Fig. 13 show Stirling engine temperature for (Plate 1 and Plate 
2) in Stirling engine container, Fig. 14 show input and output tem¬ 
perature for (Plate 1 and Plate 2), Fig. 15 show plate 1 and input 
temperature for Stirling engine container, and Fig. 16 show plate 


Stirling Engine Temperature Outputs 



—> Plate 2 Temperature Output Temperature in Plate 2 

Fig. 16. Plate 2 and output temperature for Stirling engine container during day hours on 2 March 2013 at Giza. 
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Fig. 17. Estimation of Stirling engine temperature by Matlab GUI/Interface simulation model on 2 March 2013 in Giza at 1:00 PM. 

























































638 


J.H. Shazly et al./Energy Conversion and Management 79 (2014) 626-639 


2 and output temperature for Stirling engine container during day 
hours on 2/3/2013 at Giza. 

Simulation Matlab/GUI program shows in Fig. 17 estimation of 
solar Stirling engine temperature data during day hours on 2 
March 2013 at Giza. 

6.3. Power outputs results for Stirling engine 

Figs. 18 and 19 show the results data of the solar-powered Stir¬ 
ling engine shaft power and rotating speed. Under the atmospheric 
pressure, the engine produces a maximum power of 0.5171 W at 
126 rpm and air temperature for heating container 51.52 °C. The 
engine shaft power and rotating speed also increases with increas¬ 
ing solar intensity as shown in Figs. 18 and 19. The Beale number 
and the specific power demonstrated that the performance of this 
prototype, in terms of the power output, is on par with the high 
temperature design with more simple and less costly development. 

Fig. 18 shows the Stirling engine speed during the daytime 
using Tachometer. It is shown in Fig. 18 that the maximum Stirling 


engine speed applied in Stirling engine is 126 RPM at 1:00 PM. The 
results of power data for Stirling Engine during daytime are 
showed. 

Fig. 19 shows the Stirling engine power estimation during the 
daytime using Beale formula. It is shown in Fig. 19 that the maxi¬ 
mum Stirling engine power estimation by Stirling engine is 0.5171 
watt at 1:00 PM. 

Simulation Matlab GUI/Interface model showed the shaft power 
in Fig. 20 for the solar-powered engine at different solar radiation 
intensity. 

Even though these figures are not so high, if we realize that the 
solar-powered Stirling engine is powered by a free heat source, this 
study would probably be more worthwhile as a starting point for 
the research of cheap solar-powered engines manufactured for rur¬ 
al and remote areas. The engine performance can be improved by 
increasing the precision of the engine parts and the heat source effi¬ 
ciency. The engine performance could be further increased if a bet¬ 
ter working fluid, e.g. helium or hydrogen, is used, instead of the air 
and/or by operating the engine at some degree of pressurization. 
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Fig. 18. Stirling engine speed during day hours on 2 March 2013. 
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Fig. 19. Estimation of power output from Stirling engine during day hours in Giza on 2 March 2013. 
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Fig. 20. Estimation of power output from Stirling engine by Matlab GUI/Interface simulation model in Giza at 1:00 PM. 


7. Conclusion 

This paper presents the modeling for a prototype of the solar- 
powered Stirling engine and mathematical model for the thermal 
analysis of the solar-powered low temperature Stirling engine with 
heat transfer is developed using Matlab program. The present anal¬ 
ysis provides a theoretical guidance for designing and operating of 
the solar-powered low temperature Stirling engine system, as well 
as estimating output power from the solar Stirling engine using 
Matlab Program. Also, solar radiation simulation program presents 
the implementation of a solar radiation data using MATLAB/GUI 
interface. The model is developed using basic equations of the solar 
radiation including effects of time change and different location of 
sites. 

Results from this study indicate that the engine performance in¬ 
creases with increasing solar intensity. The engine shaft power, 
speed, and heating source temperature also increase with increas¬ 
ing solar intensity. In fact, it can be said that the maximum engine 
shaft power increase with increasing heating source temperature. 
The Beale number of this engine increases with decreasing temper¬ 
ature ratio or with increasing heating source temperature. At the 
maximum solar intensity of 900 W/m 2 , or heat input of 8.55 J/s, 
the engine produced a maximum shaft power of 0.5171 W at 
126 rpm, with a heating source temperature of 64.2 °C and air tem¬ 
perature for heating container of 51.52 °C. 
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